Five-body calculation of resonance and scattering states of pentaquark system 
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Scattering problem of the uudds system, in the standard non-relativistic quark model of Isgur- 
Karl, is solved for the first time, by treating the large five-body model space including the NK scat- 
tering channel accurately with the Gaussian expansion method and the Kohn-type coupled-channel 
variational method. The two-body interaction that reproduces observed properties of ordinary 
baryons and mesons is applied to the pentaquark system with no additional adjustable parameter. 
The five-body wave function calculated has the correct asymptotic form in its the scattering-channel 
component and describes qq and qq correlations properly. The NK scattering phase shift calculated 
shows no resonance in the energy region of the reported pentaquark O + (1540), that is, at — 500 
MeV above the NK threshold (1.4 — 1.9 GeV in mass). The phase shift does show two resonances 
just above 500 MeV: a broad i + resonance with a width of V ~ 110 MeV located at ~ 520 MeV 

(~ 2.0 GeV in mass) and a sharp | resonance with V =0.12 MeV at 540 MeV. Properties of these 
states are discussed. 
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I. INTRODUCTION 

The observation of a signal of a narrow resonance at 
- 1540 MeV with S = +1 by the LEPS group 0,0, now 
called O + (1540), triggered a lot of theoretical works on 
multi-quark systems |3J, although further experimental 
reexamination is still needed to establish the state. The 
question whether the multi-quark baryon + (154O) re- 
ally exists or not is then one of current issues in hadron 
physics. In order to answer the question theoretically, 
one has to nonperturbatively evaluate the mass and the 
decay width of O + (1540), namely of the uudds resonance. 
All nonperturbative analyses made so far, however, did 
not impose any proper boundary condition to the NK 
scattering component of the pentaquark state. At the 
present stage, the non-relativistic quark model provides 
a nonperturbative framework which makes it possible to 
impose a proper boundary condition to the NK scatter- 
ing component. 

In this paper, we take the standard quark model of 
Isgur-Karl [J, y| • The Hamiltonian consists of the con- 
fining potential of harmonic oscillator type and the color- 
magnetic spin-spin interaction. As shown later, the 
Hamiltonian satisfactorily reproduces observed proper- 
ties of ordinary baryons and mesons. The same Hamil- 
tonian is applied to the pentaquark with no adjustable 
parameter. 

Resonant and non-resonant states of the uudds sys- 
tem are nonperturbatively derived with the Kohn-type 
coupled-channel variational method Q in which a proper 
boundary condition is imposed to the NK scattering 
channel and the total antisymmetrization between quarks 
is explicitly taken into account. Reliability of the method 
for scattering problem between composite particles was 
shown by one of the authors (M.K.) and actually it 
was already applied to qqq — qqq scattering ■ 

The coupled-channel variational method is accurate, 
only when the five-body dynamics in the interaction re- 
gion is solved precisely. As such a method, we use the 
Gaussian expansion method (GEM) GEM is one of 



the most reliable few-body methods proposed by two of 
the present authors (E.H. and M.K.) and their collab- 
orators. The method was successfully applied to vari- 
ous types of three- and four-body systems Q . For in- 
stance, the mass of antiproton was evaluated pre- 
cisely, i.e. with eight-digit accuracy, by comparing the 
three-body GEM calculation P3 with CERN's high- 
resolution laser spectroscopy data on highly excited 
three-body resonance states of antiprotonic helium atom 
( 4 He++ +p + e~). 

The five-body model space considered here consists of 
two parts: the asymptotic-region part describing NK 
scattering and the interaction-region part. In this paper, 
the interaction-region part is accurately described as a 
superposition of an enormous number of L 2 -type basis 
functions. As internal coordinates of the five-body sys- 
tem, we take five types of Jacobi-coordinate sets shown 
in Fig. ^ Advantages in using several rearrangement 
Jacobi-coordinate sets simultaneously are reviewed in . 
This setting can accommodate a wide (practically suffi- 
cient) model space, as shown later. Sets c = 4 and 5 
of Fig. ^ contain internal coordinates of two qq pairs, so 
these sets can treat qq correlations properly with the ba- 
sis functions of the internal coordinates. Similarly, sets 
c = 1 — 3 are convenient for handling qq and qq cor- 
relations. In general, the two-body correlation in the 
color-singlet (qq)i pair is twice as strong as that in the 
color-antitriplet ((79)3 pair due to the SU(3) color opera- 
tor. Sets 1 — 3 contain two {qq)i pairs, while sets 4 and 5 
do a (qq)i pair and a (qq)s pah- Thus, use of sets 1 — 3 is 
indispensable. In principle, four more Jacobi-coordinate 
sets are possible, but these are much less important since 
they contain only one qq or qq pair. 

The five Jacobi-coordinate sets can be classified into 
c = 1 — 3 and c = 4, 5. We call the latter two sets 
(c = 4,5) the "connected" configurations in the sense 
that they contain no color-singlet cluster; while we call 
the former sets (c = 1 — 3) the "molecular" ones as 
they are composed of color-singlet clusters. The explicit 
definitions of these configurations are shown later. Ob- 
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FIG. 1: Five sets of Jacobi coordinates among five quarks. Four u, d quarks, labeled by particle 1 — 4, are to be antisymmetrized, 
while particle 5 stands for s quark. Sets c = 4, 5 contain two qq correlations, while sets c = 1 — 3 do both qq and qq correlations. 
Sets c — 4,5 describe molecular configurations and sets c = 1 — 3 does connected ones, as shown in the text. The NK scattering 
channel is treated with c — 1. 



viously, the NK scattering component, described with 
c = 1, is molecular since it contains color-singlet qqq and 
qq clusters. 

We consider three spin-parity states, J n = h and 
J w = 2 ano - § ' W1 th a common isospin T = 0. In 
this model Hamiltonian, the J T — ^ and | states are 
degenerate in the absence of spin-orbit forces between 
quarks. (Such a spin-orbit interaction leads to mass split- 
ting between J 17 = J_ and ^ and influence on their 
phase shift Ref. [151 flfij .). Thus, our analysis is focused 
on the J 71 " = \ state. As shown later, the calculated 
phase shift exhibits no resonance when only the NK 
channel is taken and then no excitation of N and K is 
taken into account. Therefore, of importance is whether 
the five-body dynamics in the interaction region can gen- 
erate any resonance in the intermediate stage of scatter- 
ing, particularly in the energy region of 8 + (1540). Fur- 
thermore, if a resonance appears, of interest is whether 
it has the connected configuration or the molecular one. 



II. MODEL HAMILTONIAN 

Hamiltonian of the standard non-relativistic quark 
model of Isgur-Karl 0, is 



H 



Y ( mi + p- 

^\ 2m,- 



Tg + ^Conf + Vqm 



(2.1) 



where m, and are the mass and momentum of ith 
quark and Tq is the kinetic energy of the center-of-mass 
system. In what follows, u and d quarks are labeled by 
i = l—4, and s by i = 5. The confining potential Vconf 
is of harmonic oscillator type: 
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where X; is a position vector of zth quark, vq is a mass 
shift parameter, and Xf are the Gell-Mann matrices for 
color; note that A" — > — Xf for antiquark. The color- 
magnetic potential Vcm is 

-( x >-*,) 2 //3 2 tn-trj. (2.3) 
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Parameters in the Hamiltonian are fixed as follows. 
First we take standard values, m u — rrid — 330 MeV and 
m g = 500 MeV, for quark masses. And as for k and (3 
we take the same values as in Refs. 0,0, i.e. k = 455.1 
MeV- fm -2 and (3 = 0.5 fm. The remainder ^ CT and 
are so determined that the three-body calculation, done 
in the same manner as in Ref. , can reproduce tun = 
939 MeV and m A = 1232 MeV. The resultant values 
are ^/m* = -474.9 MeV and v = -428.3 MeV. The 
parameters thus determined are assumed to be universal 
for all low-lying hadrons including the pentaquark. 

One may wonder whether a linear confinement poten- 
tial should be used instead of the harmonic oscillator po- 
tential. The latter potential, however, has been used in 
many previous studies, since it is easy to handle. In our 
present analysis, rather than the practical advantage, we 
point out that as long as low lying states up to Tilo ~ 
several hundreds MeV are concerned, their properties are 
rather insensitive to the type of the potential. For this 
reason, we adopt the harmonic oscillator confining po- 
tential for the first serious calculation of the five-body 
system. (Recently, we have applied linear-type confi- 
ment and color-magnetic potential of ordinary baryons 
and mesons to the pentaquark system, and have con- 
firmed that our main conclusion of the absence of low- 
lying pentaquark resonances is unchanged in this more 
realistic calculation |l2l].) 

The present Hamiltonian is tested by static proper- 
ties of conventional baryons and mesons. Results are 
summarized in Table I which shows theoretical values of 
masses, magnetic moments and charge radii and their 
corresponding experimental values. It should be empha- 
sized that the present set of interactions satisfactorily re- 
produces observed masses of octet and decuplet baryons 
and octet mesons. We also calculated absolute strengths 
of non-leptonic weak decay of hyperon and found that 
the results reproduce the observed data owing to the qq 
correlation taken properly and are very close to those of 
the previous work 13] in which a spin-spin interaction 
different from Vcm is taken. Thus, agreement of the the- 
oretical values with the experimental values for the two- 
and three-quark systems is satisfactory to proceed to the 
five-quark system with the same Hamiltonian. 
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TABLE I: Static properties of conventional mesons and 
baryons. Squared charge radius in the last columns is denned 
by Qivf) in units of fm 2 , where Qi and r t are charge of 
the ith quark and distance of the ith quark from the center 
of mass, respectively. Magnetic moments and squared charge 
radii of states are proportional to the charge Q. 
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are the color-singlet spin ^ and spin intrinsic wave func- 
tions of ./V and K, respectively, and xl(R-i) is the wave 
function of the NK relative motion along the coordinate 
Ri (cf. c = 1 of Fig. ^) with the angular momentum 
L and the center-of-mass energy E — E t h , E t h being the 
NK threshold energy (1422 MeV in the present model). 
The reported energy of + (154O) is slightly above the 
NK threshold, it is then convenient to consider the pen- 
taquark energy E with E — E t h- 

The interaction-region part of ^j^m should be de- 
scribed with a large model space. For this purpose, we 
take five Jacobi-coordinate set, c = 1 — 5 of Fig. ^ an d 

construct the L 2 -type basis functions, ^jImoii f° r eacn 
coordinate c as follows: 



(c) 

(c) 



a = ^1234 {^(1234,5)^(1234,5) 



, sCT) (1234,5)^„ } (r c ,p c ,s c ,R c 



J" M 



} , (3.3) 



where £^ is the color-singlet wave function, and ?7o(])> 

Xs( s g<r) an< l ^L{n] are * ne isospin, spin an d spatial wave 
functions with the total isospin T = 0, the total spin 
S = h and the total orbital angular momentum L, re- 
spectively. Here we have neglected S — |, since it is 
decoupled to the N + K scattering channel with S = | . 
Here, t and s, s, a are intermediate quantum numbers of 
isospin and spin coupling, respectively, and as shown in 
the next section the symbol {n} specifies the radial de- 
pendence of the spatial wave functions. The suffix a in 



(c) 
J"M,a 



specifies a set of (t, S, s, s, a, L, {n}). 



The Hamiltonian is diagonalized within a model space 
spanned by a large number of ® jl M a , that is, ^15,000 
basis functions in actual calculations. The resulting 
discrete eigenstates are called pseudostates, when the 
eigenenergies E v satisfy E v > E t h- The pseudostates, 
{$j»m(£„);i' = 1 — i/ ma j}, are written in terms of 

*#M,a as 



III. METHOD 

We solve the five-body Schrodinger equation 
(H - E) * 7 . M = , 



(3.1) 



imposing the scattering boundary condition to the NK 
scattering component of the total wave function ^> j^m] 
here ^ j^m is classified with the angular momentum J, 
its z-component M and parity n. The NK scattering 
component of ^ j*m is expressed by 



(NK) 



— -4i234 



(E) 

{ 



^(123)^(45)] x XL (Ri)l }, (3.2 

2 2 J Jf MJ 



where the operator .A1234 antisymmetrizes particles 1 — 4 
(u, d quarks) while particle 5 is s. Here, 



and 



J*M 



£4* 

C.Q 



, Q ( r oPc 



It is possible to expand the interaction-region part of 
the total wave function in term of those eigenfunctions 
since they are considered to form a complete set for each 
J n in the finite interaction region. The total wave func- 
tion is then described as a superposition of the NK scat- 
tering component and the $jwm(E v ): 



*j«m{E) = ^™\e) 



^max 

E 

v=l 



b v {E)<S>j„ M {E v ). (3.5) 



The second term describes virtual excitations of N and K 
in c = 1 and other five-body distortions in c = 2 — 5 in the 
intermediate stage of the scattering. Unknown quantities 
Xz,(Ri) in ^jl^{E) and b v (E) are obtained by solving 
Eq. (|3.1|1 with the Kohn-type variational method |6| . In- 
clusion of high-lying pseudostates in the model space do 
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not change calculated values of the phase shift, when the 
energies E v are much larger than the NK scattering en- 
ergy. Hence, only a few tens of lowest-lying pseudostates 
contribute to the numerical calculation. In the present 
calculation, the NK scattering energy should be smaller 
than an energy of the first spatial excited state of nu- 
cleon, i.e. the Roper resonance which is located at 675 
MeV above the nucleon mass in the present model, be- 
cause we consider only the NK scattering component as 
an open channel. The condition for our calculation to be 
valid is then E — E t h < 675 MeV. For simplicity, we also 
ignore the NK* channel, since we found that inclusion 

of the channel little affects the conclusion of the present 

i + 3 
work. We also analyzed the J 71 " = | resonance of S = | 

and L = 1. It is not coupled to the NK scattering chan- 
nel with S = I but to the NK* channel. Therefore, it 
will have a small width, if it appears. However, we con- 
firmed that this is not realized at E — E t h < 675 MeV. 
These points will be discussed in the forthcoming paper. 

IV. BASIS FUNCTIONS 

Explicit definitions of the color wave functions 
^ c) (1234,5) are as follows: 

= [(123) 1 (45) 1 ] 1 , d 23 = [[(12) S (45) 1 ] S 3] 1) 
ei 3) = [(12) g [(45) 1 3] 3 ] 1 , ^ = [[(12)3(34)3] 3 5] 1 , 
# ) = [(12)3[(34)35] 3 ] 1 (4.1) 

with obvious notation; note that Q 1 ' = Q 2 ' = £j and 

(4) (5) (c) 

£i = £i due to the recombination of colors. The £j 
are connected for c = 4, 5 and molecular for other c, since 
there is no color-singlet cluster in the case of c = 4, 5. 

The isospin wave functions (1234, 5) with the to- 
tal isospin and the intermediate isospin t(= 0, 1) are 
described as 

<]) = [[(12)t3] i (45) 4 ] , „$, = [[(12) t (45) 4 ] i 3] , 

= [(12) t [(45),3] t ] Q , = [[(12) t (34) t ] 5] 0> 
^) = [(12) t [(34) t 5] t ] . (4.2) 

The spin functions v\l _ ,(1234,5) are also described as 

4U) = [[(12) S 3].(45),] X f (ssa) = [[(12).(45),] <F 3] i 
xi 3 ( U = [(12). [(45)13],] X f (ssa) = [[(12) S (34) S -] CT 5] 
4 53 sS(7) = [[(12) s [(34) s -5] ( ,] 4 . (4.3) 

Finally, the spatial wave function ipL{n}( rc ^ Pa s e> R«) 
with the total orbital angular momentum L, where L = 
for J 71 " = | and L = 1 for J" = i + , is assumed as 

= ^oo(rc) ^ooCPc) ^loo(Sc) ^ LcM (R c ). (4.4) 



Here we have set the orbital angular momenta {£ Xc } 
associated with coordinates x c = r c ,p c ,s c to be zero 
and the angular momentum L c associated with R c to 
be L; as shown below, however, this does not mean 
that the angular momentum space is small. The set 
{n} = {n r ,n p ,ns,nR} specifies the radial dependence 
of the four basis functions <fi( c \ In GEM, the functions 
</>n R L<.M(R-) are written as 

&wi.M(R) = R L ^e- (R / R "n) 2 Y LcM (R) (4.5) 

with the Gaussian ranges taken in geometric progression, 

Rn R = Rxa? 1 *- 1 (n R = 1 - ng ax ) (4.6) 

with a = (i?„max/i? 1 ) 1 /(«S ax -i). Here, n™ ax ,i?i and 
_R„max depend upon c, isospin and spin taken, but 
the explicit dependence is suppressed for simplicity 
of notation. The same procedure is taken also for 

<Coo( r c) ) <t>n,00(Pc) and 0n s ) Oo( S c) ■ 

In GEM, the model space is constructed by superpos- 
ing the Gaussian basis functions (|4.4|) over all c from 1 to 
5. This superposition is inevitable to describe few-body 
wave functions accurately, particularly when the wave 
functions have properties of strong short-range correla- 
tions and/or long-range tails; many examples are shown 
m 0. As a consequence of the superposition, further- 
more, the fact that {L Xc } = in l|4.4|) does not mean 
that the angular momentum space is small. As a sim- 
ple example, let us consider a base with {L Xcl L c } = 

0, e -(r=/f) 2 -(Pc/p) 2 ~(SJSf~(R c /Rf Xhe baS( T can bfi 

rewritten in terms of c', i.e. r c /,p c /,S c / R C '. The trans- 
formed form contains many terms with {L x >,L c i} ^ 0. 

Thus, the five-body eigenfunctions $j.jf(^) given by 
the superposition can cover an angular momentum space 
large enough to derive the phase shifts accurately. 

The antisymmetrization between quarks 1 — 4 requires 
s + 1 = even for c = 1 — 3 and s + t = s + t = even for 
c = 4, 5. For c = 1 — 3, s = 0is taken to omit the K* 
component. The diquark model proposed for J 71 " = i + in 
Ref. 01 i s included in our model space as a configuration 
of c = 4 with t — s — s = a = and L = 1, though the 
diquark is treated as a boson in the model but not in the 
present approach. 

The calculated phase shifts are converged with 
n™, n^ ax , = 5 or 6 for c = 1 - 5 and 7ig ax = 12 

for c = 1 and 6 for other c. Minimum and maximum 
ranges of the bases are, respectively, 0.2 fm and 4.0 fm 
for coordinate Ri and ~ 0.2 fm and ~ 2.0 fm for other 
coordinates. Eventually, the model space is constructed 
by ~ 15, 000 basis functions. 

It is worth noting that the connected configuration of 
c = 4, 5 has a non-negligible overlap with the molecular 
configuration of c = 1 — 3 when they are localized in a 
small space, and in general the overlap is enhanced by 
the antisymmetrization ^4i234- An extreme case is the 
(Os) 5 J v = i shell-model configuration in which the 
center-of-mass motion is excluded. The configuration is 
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equivalent to not only a molecular configuration 

A 2 34{d 1} <o) xf m f^~^~^~^} (4-7) 
but also to an connected configuration 

A xf (ln) e-£-£-^-^} , (4.8) 

where 6 = yjh/m u uj — 0.63 fm and Tilu = 

(5h 2 k/3m u )i = 300 MeV. Thus, the connected and the 
molecular configuration has 100% overlap. The (Os) 5 
wave function is obtained as an eigenstate of the approx- 
imate Hamiltonian in which m s — m u and Vcm — are 
taken and the color factor ^ a Xf A" is replaced by its 
average value —4/3. The expectation value of the full 
Hamiltonian by the (Os) 5 J 71 — \ configuration is 503 
MeV above the NK threshold. The (Os) 5 configuration 
is not an eigenstate of the full Hamiltonian and has a 
large overlap with the NK scattering configuration, so it 
is changed into non- resonant continuum states, as shown 
later, when the Schrodinger equation is solved with the 
NK scattering component. 

V. RESULTS 

Firstly, we calculate the NK elastic scattering phase 
shifts for J" = i and ^ + , by omitting the pseudostate 
terms in <|3.5I) . In this calculation, N and K behave as 
inert particles. The resulting phase shifts are shown in 
Fig. |21 as dash-dotted lines. There appears no resonance 
in the entire energy region. 

Next, we do the full-fledged calculation including the 
pseudostate terms in <|3.5[1 . Calculated phase shifts are 
shown as solid lines. Up to the energy E ~ 300 MeV, the 
calculated phase shifts of J 7 * = 1/2 " agree qualitatively 
well with the experimental data [la . Ila| as well as with 
the previous quark model estimations (l7|. 

No resonance is seen in the energy region — 500 MeV 
above the NK threshold (1.4— 1.9 GeV in mass), that is, 
in the reported energy region of O + (1540). One does see 

two resonances around 530 MeV; one is a sharp J n = | 
resonance with a width of Y = 0.12 MeV located at 
E — E t h — 540 MeV and the other is a broad | res- 
onance at ~ 520 MeV with Y ~ 110 MeV. It should be 
noted here that our result on the absence of the low-lying 
pentaquark resonance is consistent with some recent lat- 
tice Q CD results in Refs. jHlilEl, 

in which, in fact, no 
pentaquark resonance is observed below about 300 MeV 
with respect to the N + K threshold. 

We investigate properties of the sharp J n = i res- 
onance by using Eq. I|3.5|) in its second term. In the 
interaction region which we are interested in, the am- 
plitude of the second (pseudostate) term is much larger 
than that of the first (scattering-channel) one. Thus, the 
sharp resonance is not a molecular state of N and K in 
their ground states. Figure shows the one-body den- 
sity of the resonance state in its u,d quark component 
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FIG. 2: Calculated phase shifts for a) = | and b) 

J n = | + states. The solid lines are given by the full-fledged 
calculation, while the dash-dotted lines are by the approx- 
imate calculation with the elastic NK channel alone (see 
Encr gies are measured from the NK threshold. The 
arrow indicates the energy of + (154O) in E — E t ^. 



and s quark one. One can see from the density that a 
s quark is localized near the center and surrounded by 
u, d quarks which are totally antisymmetrized but spa- 
tially symmetrized. The density of u, d quarks is found 
to be very broad compared with that of the (Os) 5 state. 
This indicates that u, d quarks are not in (0s) shell but 
in higher shells. Actually, the overlap of the second term 
of Eq. (|3.5(l with the (Os) 5 configuration is only 2 %. 
Furthermore, the one-body r.m.s. radius measured from 
the center of mass is 1.10 fm for u, d quarks and 0.72 fm 
for s quark, while the corresponding radius is 0.69 fm 
for the (Os) 5 configuration. Thus, the structure of the 
sharp resonance is quite far from the (Os) 5 configuration, 
although the expectation energy of the (Os) 5 configura- 
tion is rather close to the resonance energy. As for the 
broad J n = \ + resonance, the one-body density cannot 








FIG. 3: One-body density p(r) of the J 71 = \ resonance at 
E — E t h = 540 MeV as a function of r that denotes the dis- 
tance of a quark with respect to the center of mass. The den- 
sity of u, d quark is given by the solid line while that of s quark 
is by the dashed line. p(r) is normalized as J °° p(r)r 2 dr = 1. 



be estimated meaningfully, since the second term are not 
dominant in Eq. (|3.5|l . 

In general neither pure connected nor pure molecu- 
lar resonance exists, since any connected and molecular 
configurations are non-orthogonal to each other due to 
the antisymmetrization. Nevertheless, we can determine 
whether the two resonances are connected or molecular 
in their main components. The phase shifts shown above 
are not changed significantly when the connected com- 
ponents (c=4, 5) are omitted from the basis l|3.3[) : for ex- 
ample, the resonance energy goes up only by 6 MeV for 



the J" 



state and by 15 MeV for the J*" = 



state. 



In contrast, when the molecular components (c = 1 — 3) 
are omitted, the J 71 = i resonance disappears from 

the entire energy region and the J™ = i + resonance is 
shifted upward by 130 MeV. This result means that the 
two resonances are molecular in their dominant compo- 
nents. Thus, as for both the resonances the five-body 
configurations in the interaction region, represented by 
the second term of (|3.5|l . are mainly molecular, that is, 
they are composed of color-singlet but spatially distorted 
(excited) (qqq)i and (qq)i clusters. The reason why the 
five-body system is mainly molecular at energies con- 
cerned in this work is as follows. Since the (qq)± cor- 
relation is twice as attractive as the (qq)s one, in general 
the molecular configuration with one qq pair and one qq 
pair (cf. c = 1 — 3 in Fig.^l obtain an energy gain larger 
than the connected configuration with two qq pairs (cf. 
c = 4, 5) does. Thus, the molecular five-body config- 
uration in which the qq and qq correlations work most 
effectively appears as low-lying states of the pentaquark 
system. For the two resonance states {qqq)\ and (qq)i 
clusters are spatially distorted to a large extent, while 
for non-resonance states the distortion is rather weak. 



Exotic resonances might appear at energies much higher 
than the two molecular resonances do, but it is out of 
scope of this paper since some inelastic scattering chan- 
nels are opened there. 



a) J= 



600- 



% 400- 



LU 

I 
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... : 
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+ 



1.6 1.8 



FIG. 4: The stabilization plots of the eigenenergies E v of the 
pseudostates against the scaling factor a for a) J n = i and 



b) r 



Here, the Gaussian ranges Rn R associated with 



the coordinate Ri is scaled as R nR — ► aR nR with 0.7 < a < 
1.8; The arrow indicates the energy of O + (1540) in E v — Eth.- 

All the nonrelativistic quark-model calculations done 
so far for the pentaquark system use the bound-state ap- 
proximation in which the NK scattering component does 
not exist asymptotically. Therefore, it is of interest to 
compare the scattering phase shifts with the eigenener- 
gies (E v ) of the pseudostates as a result of the bound- 
state approximation. In general, in the case of a sharp 
resonance, there is an eigenenergy E v of the Hamiltonian 
that is close to the resonance energy, and the energy E v 
is stable with respect to scaling the ranges of the ba- 
sis functions associated with Ri the coordinate of the 
NK scattering, since the asymptotic amplitude of the 
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sharp resonance is relatively much smaller than the in- 
ternal amplitude. In contrast, other pseudostates have 
eigenenergies E v which decrease monotonously with the 
increasing the ranges, because so does the kinetic energy 
associated with R4. 

These characteristics are really seen in Fig. 0] that 
shows stabilization plots of E v with respect to the scal- 
ing R nR — > aR nR ; the parameter a corresponds to the 
volume size in the lattice calculation. This plotting to in- 
vestigate resonances is called the real-scaling (stabiliza- 
tion) method 18]. In Fig. Hi, for J n = |~, all pseu- 
dostates except one are unstable in the sense that they 
decrease toward the NK threshold. Thus, those pseu- 
dostates are regarded as a discrete representation of the 
non-resonant continuum spectrum. The stable (horizon- 
tal) line at E v — 540 MeV precisely corresponds to the 
sharp resonance at 540 MeV in Fig. [SJi. The horizon- 
tal line and the nearby unstable lines do not cross each 
other because of repulsive forces working between them, 
though it is not precisely plotted. This is called the avoid- 
ing crossing. One can roughly estimate the width of the 
resonance from the behavior of the avoiding crossing with 
the aid of Eq. (4) of Ref. ^| . In the present case, the es- 
timated width is of order 0.1 MeV and consistent with the 
value obtained from the phase shift. On the other hand, 
the ,V = \ + resonance at ~ 520 MeV with Y ~ 110 MeV 
in Fig. \Bp is too broad to be identified as a resonance 
from the stabilization plot of Fig. 0Jd, though one sees a 
tendency of plateaus around 500 — 600 MeV. Therefore, 
without the scattering calculation it is difficult to dis- 
criminate a broad resonance from non-resonance states. 

In summary, we solved the five-quark scattering prob- 
lem by applying GEM |8| and the Kohn-type variational 
method p to the large model space including the NK 
scattering component. We adopted the standard non- 
relativistic quark model of Isgur-Karl which satisfactorily 
reproduces experimental values of the two- and three- 
quark systems. The resultant NK scattering phase shift 
showed no resonance in the reported energy region of 
G + (1540); this is the most important result of the present 
paper. At energies much higher (by ~ 400 MeV) than the 

O + (1540) energy, we did find a broad J v = \ resonance 
with r ~ 110 MeV at ~ 520 MeV above the NK thresh- 
old and a sharp J* — | resonance with T = 0.12 MeV 
at 540 MeV. In the present model Hamiltonian, since the 
J" = i + and | + states are degenerate to each other, 
there exists a broad J 7r = | resonance with the same 



energy and width as the J 71 " = \ resonance. 

We have done the same calculation for other Hamilto- 
nian proposed in [T^l ] that also reproduces the observed 
properties of the ordinary hadrons and mesons with the 
same quality of agreement as in Section II. The result 
was qualitatively the same as in this paper; in partic- 
ular resonances are absent in the low energy region up 
to 500 MeV from the NK threshold. The locations and 
the widths of the resonances at higher energies, however, 
depend on the details of the model hamiltonian. 

The resonance states are mainly composed of color- 
singlet (qqq)i and (qq)i clusters which are distorted (ex- 
cited) to a large extent. Thus, these are molecular res- 
onances. The {qq)i correlation is twice as attractive as 
the (qq)s one, so low-lying states, no matter whether the 
widths are small or not, are dominated by molecular con- 
figurations. Thus, even if a resonance with a small width 
is measured, the fact that the width is small does not 
necessarily mean that the resonance is connected. The 
sharp J 71 " = i resonance has a quite different structure 
from the (Os) 5 configuration. In the resonance, s quark 
is located near the center and surrounded by u, d quarks 
which are spatially symmetrized. 

Finally, we tested the reliability of the bound-state ap- 
proximation with the real scaling method. Essentially the 
same approximation is used in the lattice calculation. We 
found, by comparing Fig. 4 with Fig. 2, that the approx- 
imation surely works for a sharp resonance with a width 
of order 0.1 MeV, but not for a broad resonance with 
a width of order 100 MeV and that most of the states 
obtained by the approximation melt into non-resonant 
continuum states when the NK scattering channel is in- 
cluded accurately. It should be noted that our model 
calculation with the NK scattering channel clarifies the 
mechanism of the disappearance of low-lying pentaquark 
resonace, which cannot be shown even by lattice QCD. 
Further analyses including NK* scattering channel will 
be reported in a forthcoming paper, together with results 
for other (T, J 77 ) states. 
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